Charge redistribution at the antiferromagnetic phase transition in SrFeAsF compound 
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The relationship between spin, electron, and crystal structure has been one of the foremost issues in 
understanding the superconducting mechanism since the discovery of iron-based high temperature super- 
conductors. Here, we report M(3ssbauer and first-principles calculations studies of the parent compound 
SrFeAsF with the largest temperature gap (~50 K) between the structural and antiferromagnetic (AFM) 
transitions. Our results reveal that the structural transition has little eff^ect on the electronic structure of the 
compound SrFeAsF while the development of the AFM order induces a redistribution of the charges near 
the Fermi level. 

PACS numbers: 76.80.-Hy, 74.10. -nv, 71.15.Mb 



INTRODUCTION 



EXPERIMENTS 



The discovery of superconductivity (SC) in iron-based 
superconductors, like the cuprates, promises to be an im- 
portant milestone in condensed matter physics 1 1-3]. Prox- 
imity to a magnetically ordered state has been suggested to 
be a crucial factor for the observed high-temperature (high- 
Tc) SC 10] ■ Another interesting feature is that the anti- 
ferromagnetic (AFM) transition is always accompanied by 
a structural transition either before the AFM transition or 
at the same temperature (see [5] and references therein). 
In the 122-type AeFe2As2 (Ae=Ba, Sr, Ca, Eu) parent 
compounds and the 11-type Fei+^Te systems, the struc- 
tural and AFM transitions occur simultaneously. For the 
rare-earth based 1 1 1 1 -type ReFeAsO (Re=La, Ce, Pr, Nd, 
etc.) materials and the Ill-type AFeAs (A=Li, Na) com- 
pounds, the structural transition precedes the AFM transi- 
tion by about 10~20K. Till now, the relationship between 
the structural and magnetic phase transition has been an 
open question though it has been suggested that the struc- 
tural transition may be driven magnetically [ 6] . But there is 
angle-resolved photoemission spectroscopy (ARPES) evi- 
dence [7] that structural and AFM transitions may share 
the same origin and could both be driven by the electronic 
structure reconstruction in NaFeAs. Actually, the nature of 
the AFM order itself is still under debate, with some favor- 
] ing a localized picture and others supporting an itinerant 
nature. 

I To understand the interplay between structural, mag- 
netic, and electronic properties, a microscopic investiga- 
tion by Mossbauer spectroscopy may be very help- 
ful. In the present work, we report detailed Mossbauer 
spectroscopy measurements and first-principles calcula- 
tions of SrFeAsF. The SrFeAsF parent compound, which 
has the largest reported temperature interval (~50K) be- 
tween the two phase transitions |14], offers an ideal op- 
portunity to address what happens between the structural 
and AFM transitions. Our results show strong evidence of 
charge redistribution at the AFM phase transition in the Sr- 
FeAsF compound. 



SrFeAsF was synthesized using a two-step solid-state re- 
action method, which is similar to that described earlier 
ifisll . First, SrAs/FeAs was prepared by reacting Sr/Fe 
flakes/powders and As powders sealed in an evacuated 
quartz tube at 773 K for 10 h and then at lOOOK for 20 h. 
The resultant precursors, SrAs, FeAs, SrF2 and Fe pow- 
ders, were mixed and thoroughly ground at a molar ratio of 
1:1:1:1. All the weighing and mixing procedures were per- 
formed in a glove box with protective argon atmosphere. 
Then the mixture was pressed into a pellet, sealed in an 
evacuated quartz tube and heated up to 1273 K at a rate of 
60 K/h, kept at this temperature for 40 h, and cooled down 
to room temperature. The product was homogenized in an 
agate mortar, pressed into a pellet, and sintered at 1273 K 
for 40 h again to obtain the final sample. 

Phase purity was checked by x-ray powder diffraction 
(XRPD) on a Philips X'pert diffractometer with Cu Kq, 
radiation. Rietveld refinement of the SrFeAsF was per- 
formed with the GSAS package liT6ll to obtain the lat- 
tice parameters. Transmission Mossbauer spectra (MS) at 
temperatures between 16 and 290 K were recorded using 
a conventional constant acceleration spectrometer with a 
y-ray source of 25 mCi ^^Co in a palladium matrix mov- 
ing at room temperature. The absorber was kept static in 
a temperature-controllable cryostat filled with helium gas. 
The drive velocity is calibrated with sodium nitroprusside 
(SNP) at room temperature and all the isomer shifts quoted 
in this work are relative to that of the a-Fe. 



CALCULATION METHODS 

We use the full potential linearized augmented plane 
wave (FP-LAPW) iITtIi method as embodied in the 
WIEN2K ilil 11911 code in a scalar relativistic version 
without spin-orbit coupling. According to the FP-LAPW 
method, the simulation cell is divided into spherical atomic 
regions with radii R,„, (2.30, 2.41, 2.13, and 2.30 a.u. for 
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Sr, Fe, As, and F atoms, respectively.) and an interstitial 
region. The wave functions are described by atomic like 
functions inside the atomic spheres, while in the intersti- 
tial region plane waves are used. The maximum angu- 
lar momenta / for the expansion of the wave functions in 
the spherical harmonics inside the spheres is confined to 
Imax = 10. The cutoff parameter Rmt^max for limiting the 
number of the plane waves is set as 7.0. The magnitude of 
the largest vector in the Fourier expansion of the electron 
density (Gmajr) is set to 14.0. States lying more than 6Ry 
below the Fermi level are treated as the core states. The 
Brillouin zone integration is done with a modified tetrahe- 
dron method using 648 special ^-points in the irreducible 
wedge. The generalized gradient approximation (GGA) 
suggested by Perdew, Burke and Emzerhof (PEE GGA) 
ll20ll is employed for the exchange-correlation effects. 

Once the electron (spin-) densities are calculated self- 
consistently with high accuracy, hyperfine parameters can 
be deduced directly ll2ll l22ll . The electric field gra- 
dient (EFG) tensor can be obtained from an integral 
over the nonspherical electron density. The most in- 
teresting principal component V„ is given by V,, = 

2^-0 Iq" £ 0^2 (cos 6) f p(r), where P2 is the 
second-order Legendre polynomial and p(r) is the elec- 
tron density. In the LAPW method, p(r) is expressed 
by lattice harmonics, i.e. under an LM form, plm(j) = 
ZEE ^/,„(r)i?„„,(r)G^,"'«'. It turns out, only the 

E;,f;<Ef l,m l\m' 

component p2o is important to V„. And when L = 2 and 
M = 0, only I = V = I and / = /' = 2 (and to a minor 
extent I = 0, F = 2 and / = 1 , /' = 3) give non-zero Gaunt 
numbers, which correspond to p-p and d-d contributions to 
the EFG, respectively. 



RESULTS AND DISCUSSION 

Figure [1] shows the XRPD pattern of the sample, from 
which one can see that all the main peaks can be indexed 
to the SrFeAsF phase with a ZrCuSiAs-type structure. The 
lattice constants were determined to be a=4.002A and 
c=8.970 A, which is in good agreement with previously re- 
ported values [T?]. Only a small amount of the SrF2 im- 
purity phase was detected, which is determined by the Ri- 
etveld refinement to be less than 4%. As the impurity phase 
(SrF2) does not contain iron, the analyses of the MS can not 
be affected and this is shown in what follows. 

M5ssbauer spectra recorded at indicated temperatures 
together with transmission integral fits are shown in Fig. 
Ill Clearly, the spectra change from the paramagnetic (PM) 
state (singlet or unresolved doublet) to the magnetically 
split state (sextet) upon lowering of the temperature. The 
Hamiltonian for the hyperfine coupling in the principal axis 
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FIG. 1 . (Color online) X-ray powder diffraction pattern of Sr- 
FeAsF (♦) and Rietveld fit (red solid line). Bars at bottom show 
the calculated Bragg diffraction positions of SrFeAsF (upper row) 
and a small amount of SrF2 (lower row) impurity phase (less than 
4%). The olive green line is the difference between the observed 
and calculated patterns. 



system of the EFG can be expressed as 12311 



M 



(1) 



4/(27-1) 

-g/u^BKIx cos (p + ly sin (J) sin 6 + 1^ cos 6\, 

where 7 and Ix, ly, Iz refer to the nuclear spin operator and 
operators of the nuclear spin projections onto the princi- 
pal axes and Q denotes the quadrupole moment of the nu- 
cleus (for excited ^^Fe Q = -i-0.17(l)b H). and 9 
are the azimuthal and polar angles of the hyperfine mag- 
netic field in the EFG coordinate system, respectively. 
V = (^xx - Vyy)/Vzz represents the asymmetry parameter 
of the EFG at the nucleus. In the present case, it is reason- 
able to assume that the EFG has axial symmetry around 
the c-axis ifToll . Therefore, above the AFM transition 
temperature, one has a doublet with quadrupole splitting 
AEq = j\eQVzz\- Below the AFM transition temperature, 
due to the smallness of the quadrupole interaction in com- 
parison with the magnetic dipole interaction, the assump- 
tion of first-order perturbation can be made. Hence, eigen- 
values of the magnetic Hamiltonian are perturbed (shifted) 
by the term e = {-l f"'^^^'^ieQVzz/S) ■ (Scos^B- 1) [23]. 
Considering the established magnetic structure, 9 = 90°, 
the sign of Vzz and the nuclear quadrupole coupling con- 
stant (NQCC), \eQVzz.\, can be obtained subsequently. 

As can be seen, the spectra can be well fitted by only one 
double t/(sextet) above/(below) the transition region, indi- 
cating that the local environment of the Fe ion is unique. 
This means that there is no Fe containing impurity phase 
in our sample, coincident with the XRPD results. At 
room temperature the best fit of the spectra gives an iso- 
mer shift of S =0.44 mm/s and a quadrupole splitting of 
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\AEq\ -0A2 mm/s, which compares well with reported val- 
ues for SrFeAsF fl^ and other 1111-type iron-pnictides 
ifioi [Till . At temperatures in the transition region (100 K- 
125 K) the MS are fitted with the superposition of a doublet 
and a sextet iflollTlll as usually done in the Mossbauer stud- 
ies of the iron-based parent compounds. 
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FIG. 2. (Color online) Mossbauer spectra taken at indicated tem- 
peratures of the SrFeAsF sample together with transmission inte- 
gral fits (solid line). 

The temperature dependence of the hyperiine field, 
Bf,f(T), is shown in Fig. [3] Fitting Bijf(T) data to the usual 
function fi/,^(r) = 5/,/(0)[l -(r/rA,)"f LI4J (reddashdot- 
ted line) leads to 5,,/(0)=4.84T, 7^=125 K, a = 2.8 and 
yS = 0.23. These values agree very well with those derived 
from fitting the upper precession frequency of the jiSR data 
lfl4ll and prove the two-dimensional magnetism ll26ll in the 
SrFeAsF compound. To find out the order of the magnetic 
phase transition, we fit the Bi,f(T) data according to the 
Bean-Rodbell model (blue solid line) IHIH for 7 = 1/2, 
T/Tq = (cr/ tanh"' cr)(l -i- where Tq is the tran- 

sition temperature if the lattice is not compressible, cr is 



the sublattice magnetization, and ^ is a fitting parameter 
(^ < 1 for a second order transition, ^ > 1 for a first order 
transition and for ^ = the equation reduces to a Bril- 
louin function). The fitted values for f , To and Bhf{T) are 
found to be 0.8, 125 K and 4.82 T, respectively. Obviously, 
f is a little smaller than 1 , which is different from that ob- 
tained for CaFe2As2 (^ = 1.35), indicating that the AFM 
transition for SrFeAsF is more second-order-like than that 
for CaFe2As2, which is in agreement with previous reports 

IF 
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FIG. 3. (Color online) Temperature dependence of the hyper- 
fine field, BhfiT), extracted from transmission integral fits of the 
Mossbauer spectra. The red dashdotted line and the blue solid 
line are fits to the data with the power law and the Bean-Rodbell 
model, respectively (see text). 

Temperature dependence of the isomer shift, 6{T), de- 
termined from the fits of the MS is shown in Fig. [4] (a). 
The two vertical dashdotted lines at 125 and 175 K lfl4ll 
indicate the AFM and structural transition temperatures, 
respectively. As expected, 5{T) increases gradually with 
decreasing temperature due to the second order Doppler 
effect. And in terms of the Debye model, the Debye tem- 
perature is found to be 0£)=337K, which is in reasonable 
agreement with those determined from specific-heat mea- 
surement, 0o=339K 1251], and neutron diffraction experi- 
ment, 0£)=347 K i28n . In Fig. |4](b) we present the temper- 
ature dependence of the NQCC, \eQV--\, which is obtained 
as 2AEq and 8|e| for PM and AFM states, respectively. 

Obviously, no abrupt changes of 5(7) and NQCC dur- 
ing the structural transition around ~175 K are found, im- 
plying a rather continuous structural transition. In terms of 
group theory, the space group of the low-temperature phase 
(Cmma) is a subgroup of that of the high-temperature phase 
(P4/nmm) and therefore a relatively continuous structural 
transition is usually the case |25]. Considering the small 
change in volume after the transition I25L i28ll . it is likely 
that the charge density at the iron nucleus and the spa- 
tial distribution symmetry of the charges exhibit negligible 
changes. Therefore, one may not expect distinct changes 
in 6(T) and NQCC near the structural transition tempera- 
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FIG. 4. (Color online) Temperature dependence of the isomer 
shift S together with theoretical fits using the Debye model, (a), 
and the nuclear quadrupole coupling constant, (b), of the Sr- 
FeAsF compound. The blue-diamonds (purple-circles) are de- 
rived from the AFM (PM) phase. Vertical dash-dotted lines at 
125 and 175 K [14,1 indicate the AFM and structural transition 
temperatures, respectively. 



ture range. Interestingly, upon cooling through the AFM 
transition, the NQCC jumps to a much smaller value for 
the AFM phase while the temperature dependence of S(T) 
still follows the Debye model without any evident jump. 
The isomer shift changes rather smoothly with decreas- 
ing temperature through the AFM transition was also ob- 
served in the GdFeAsO compound 111 111 in contrast to that 
for the 122-t ype CaFe2As2 material where a small jump 
was detected il2ll . This could be understood as due to some 
charge redistribution effect, which changes the EFG at the 
iron nucleus while has little effect on the absolute charge 
density at the iron nucleus. In fact, charge redistribution 
effect have also been observed in 122-type AeFe2As2 ll29ll 
and Ill-type NaFeAs liTJ systems and high-Tc cuprates 
materials 

To get a better understanding of the charge redis- 
tribution effect, we calculate the EFG and the elec- 
tronic properties of the SrFeAsF compound. For tetrag- 
onal (P4/nmm) SrFeAsF, the experimental lattice pa- 
rameters and internal atomic positions (a=b=3.9930 A, 
c=8.9546A, Z5,=0.1583, Za5=0.3485) are adopted iH. 
For orthorhombic (Cmma) SrFeAsF, the experimental 
value (a=5.6155 A, b=5.6602 A, c=8.9173 A, Z5, =0.1584, 



ZAi=0.3475) are employed UJ]. The AFM phase of Sr- 
FeAsF is also set up according to experimental reports, that 
is, the spin direction is opposite between two adjacent FeAs 
layers and opposite along the longer Fe-Fe direction within 
the FeAs layers [28]. 

First, we have calculated the EFG principal component 
Vj., for nonmagnetic -P4/nmm (NM-P4/nmm), NM-Cmma 
and AFM-Cmma phases of the SrFeAsF parent compound. 
The experimental and numerical results are shown in Ta- 
ble I. The sign of V,, for NM-P4/nmm and NM-Cmma 
phases can not be determined experimentally and the ab- 
solute values are given. For the magnetically split AFM 
state, the sign of V-,-, can be established directly in this case 
to be negative, which agrees well with first-principles cal- 
culations. This is different with other iron-based supercon- 
ductors, where the sign of Vj.j. is obtained to be positive 
1 13- 121. The reason for this difference certainly deserves 
future investigations. As can be seen, if we consider the ab- 
solute values of V-,, the tendency of the calculated results 
match reasonably well with the experimental results. The 
\V,,\ values are almost the same for NM-P4/nmm phase and 
the NM-Cmma phase. While the |V--| value for the AFM- 
Cmma phase is significantly smaller than the value for the 
other two phases. Generally, the EFG is sensitive to struc- 
tural changes and is taken as a good probe for structural 
transitions. Our observations, however, indicate that the 
structural transition from the P4/nmm phase to the Cmma 
phase in SrFeAsF has no obvious effect on the EFG, but the 
introduction of long-range magnetic order unambiguously 
reduces the |V,,| value of SrFeAsF. 



TABLE I. Principal component of the EFG, V^^, for NM- 
P4/nmm, NM-Cmma and AFM-Cmma phases of the SrFeAsF 
compound. The experimental values for NM-P4/nmm and NM- 
Cmma phases are given as |V;;| (see text). 

Vzz (10^' V/m^) NM-P4/nmm NM-Cmma AFM-Cmma 



Calculation 
Experimentation 



-0.350 
0.65(3) 



-0.330 
0.64(4) 



-0.197 
-0.27(1) 



In order to further understand the V,,, the evolutions of 
V^"'', Vff^ and V^^ with the state of energy are investi- 
gated minutely and the results are shown in Fig. |5](a). For 
all three structures, V^ '', V^~'' and V„ are nearly zero un- 
til around -5 eV starting from the lower energy side, after- 
which Vt~'' and V^~'^ vibrate strongly until the Fermi en- 
ergy is reached. Obviously, the V^z of all systems are domi- 
nated by the valence electrons near the Fermi surface. And 
the similarity in the evolution of V-,-, with energy is quite 
distinct for the two NM systems, which explains the ex- 
perimental observation quite well. It is found that, though 
both p-p interaction and d-d interaction contribute to the 
Vzz, the contribution from V^~^ dominates. 

According to the above definition, the EFG is decided 
by the non-spherical electron-density p2o(r), as well as by 
the radial dependence 1 /P of the anisotropic charge den- 
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FIG. 5. (Color online) Evolutions of V^r'', V^'', and V,,, (a) and, 
evolutions of Awy, and An^ (b) with the state of energy (see text). 



sity. We calculated the nonspherical numbers of p and d 
charges inside the Rmt of the Fe atom, denoted as An,, 
and Arid, respectively. Atip and Arid quantify the devi- 
ations of the charge distribution from spherical symme- 
try and are calculated as: Arip = (rip^ + tip^}/2 - rip^ and 
Arid = {fid,, + "c/,2_,,2) - ("rf,,- + nd,Sl2 - rid^^ HI]. The 
evolutions of Arip and Arid with the state of energy are 
shown in Fig. [5] (b). Obviously, the fluctuation of Arid 
is much more pronounced than that of Arip, indicating 
stronger anisotropic spatial distribution of ¥e-3d electrons 
compared with that of the Fe-3p electrons. This can be 
explained by the strong hybridization between Fe-3p and 
As-45 orbits, which reduces the spatial anisotropy of the 
Fe-3p charges. The decrease of Arid in the AFM phase can 
be understood from the calculated density of states (DOS) 
of the iron atom for the three structures. As shown in Fig. 
m the partial DOS of Fe-3ii sub-orbitals for the two NM 
phases almost has the same graphics features in the energy 
range from -5.0 eV to the Fermi level. Though there ex- 
ists some diff'erences above the Fermi level, the discrepan- 
cies do not affect the EFG. The diff'erences in the DOS of 
the AFM state are quite distinct. Obviously, the dxy, dxz. 



and dxz orbitals control the symmetry of the Fermi surface 
in NM phases, yet the DOS is much reduced in the AFM 
phase, resulting in a much smaller Fermi surface |33]. Be- 
sides, the gravity centers of all five d orbitals are closer to 
each other in energy, yielding more isotropic spatial distri- 
bution of d electrons. This will reduce the anisotropy of 
charge density of d electrons, hence resulting in a smaller 
EFG, proving the above observation. 




FIG. 6. (Color online) Density of states of the iron atom, (a) 
NM-P4/nmm, (b) NM-Cmma and (c) AFM-Cmma as described 
in the text. 

The fact that V^^ of the three systems is dominated by d-d 
interactions is the key point to understand the sudden drop 
in the EFG of the MS observations. It is well known that 
the magnetic moment of the AFM-Cmma phase mainly 
comes from the spin polarization of ¥e,-3d electrons. Due 
to the exchange splitting of the 3J-band, the chemical po- 
tential of 3(i-electrons would decrease in a magnetic phase. 
Consequently, there should be a redistribution of electrons 
between sp- and J-bands. The charge redistribution due 
to spin polarization of Fe-3J electrons may have at least 
two effects on the electronic structure of the AFM phase. 
First, the charge will transfer from Sr-F layers to Fe-As 
layers. Second, the 3d electrons will redistribute among 
the 'id suborbitals. 

The calculated electron density differences of the NM- 
Cmma phase and the AFM-Cmma phase prove our sugges- 
tions. Figure |7] (a) displays the electron density difference 
Ap for the two phases in the plane defined by the Fe atom 
and the two nearest As neighbors. The direction of the 
nearest Fe atoms is along the crystalline fe-axis. The pro- 
file of the density difference along the Z^-axis is also shown 
in Fig. |2](c), where a positive number denotes increased 
density after magnetization and a negative number denotes 
the decreased density. As can be seen, the most distinct 
increase in charge density happens at the iron sites, indi- 
cating the electron charges tend to accumulate at the iron 
sites after the emergence of magnetic ordering. We also 
calculated the electron density difference of the yz-plane 
with As and Sr atoms. As shown in Figs. |7](b) and (d). 
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the electrons distributed in the region between two adja- 
cent Fe-As layers decrease after magnetic ordering, which 
means the electrons tend to flow from Sr-F layers to Fe-As 
layers accompanying the emergence of magnetic ordering. 



by the EFG. These findings may also hold for other iron- 
pnictides and may shed light on understanding the fascinat- 
ing physics of these materials. 
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FIG, 7, (Color online) The difference in electron density distri- 
bution between the NM-Cmma phase and the AFM-Cmma phase 
of the SrFeAsF compound in the Fe-As plane (a) and the Sr-As 
plane (b) and the profiles of the density difference along the [010] 
direction (c) and the [001] direction (d) of the crystal. 

Moreover, it is found that the density difference distri- 
bution around the iron site clearly has spatial characters of 
3d orbitals, suggesting the redistribution of the electrons 
among the 3d orbitals. From the calculated occupation 
number of each orbital, it is found that electrons move from 



J,2, d 



^j. and dy^ orbitals to dxy and dx 



orbitals when the 



system goes into the magnetic state. Indeed, the electron 
density projected on the Fe-As plane increases along the 
b-axis after magnetization, as shown in Fig. |7](a). While 
the electron density along the direction of the As-Fe bond 
decreases. This observation highly suggests that magne- 
tization will enhance the localization of the 3d electrons, 
hence the correlation among 3d electrons. 

Actually, magnetization induced charge redistribution is 
common in terms of thermodynamics. But it is not easy 
to observe experimentally in iron based superconductors 
10, 11^ due to the interference of the structural phase transi- 
tion, which happens close to the magnetic phase transition 
and may have magnetic origin [i^. Luckily, in the SrFeAsF 
compound, the large temperature interval between the 
structural and magnetic phase transitions (~50 K) makes it 
possible to exclude the disturbance in electronic redistribu- 
tion due to the structural phase transition. For the SrFeAsF 
compound, Fe-3d electrons contribute to the EFG as well 
as to the magnetism. And the EFG is sensitive to the spatial 
distribution of charges near the iron nucleus; accordingly, 
the charge redistribution around the iron site, due to the 
appearance of long-range magnetic order, can be detected 



CONCLUDING REMARKS 

In summary, we have studied the SrFeAsF parent com- 
pound by Mossbauer spectroscopy and first-principles cal- 
culations. Our results about the structural and AFM tran- 
sitions are in good agreement with previous reports. The 
temperature dependence of the isomer shift has no anoma- 
lies during the two phase transitions, indicating no abrupt 
changes of the absolute charge density at the iron nucleus. 
Additionally, the structural transition has little effect on the 
main component of the EFG, V~~, while during the AFM 
phase transition V„ jumps to a much smaller value in the 
low temperature AFM state. Our first-principles calcula- 
tions reveal that the contribution to the EFG mainly comes 
from the electrons close to the Fermi level. These obser- 
vations could be well understood by the magnetically in- 
duced redistribution of the charges near the Fermi level, 
which changes the electron distribution symmetry (hence 
the EFG) while has little effect on the absolute charge den- 
sity at the iron nucleus during the AFM transition. 

This work was supported by the National Natural Sci- 
ence Foundation of China under Grants No. 10975066. 
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